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ABSTRACT 

A theoretical Investigation Is made of the evolution of a vapor bubble 
for a submerged Journal bearing under dynamically loaJed conditions. The 
solution to the Reynolds equation Is determined numerically using a control 
volume method (Elrod algorithm). This method conserves mass throughout the 
computational domain Including the liquid-vapor Interface which may or may not 
be In motion relative to the minimum film line. An ADI (Alternating Direction 
Implicit) method Is used to effect the time march. Excellent agreement was 
found with the experimental work of Jakobsson and Floberg for stationary 
cavitation. Predictions of bubble life for nonstationary cavitation compare 
reasonably well with that measured by Jacobson and Hamrock using high-speed 
photography. A comparison study was performed to determine some of the 
consequences of applying a nonconservative theory to a dynamic problem. A 
complete dynamic cycle of a Journal whirling In a circular path was chosen for 
the basis of comparison. Significant differences were observed In the load 
components near the end of the cycle. In each case, onset of cavitation was 
observed followed by bubble growth and subsequent collapse. More complete 
details of this phenomena are Illustrated with the use of perspective graphic 
plots depicting the associated pressure distribution and region of cavitation 
with position and motion of the Journal within the housing. 


NOMENCLATURE 

0 shaft diameter, m 

ti eccentricity, m 

Fr radial load component, N 

Ftp tangential load component, N 

F m friction force due to shear stress, N 

g switch 'unction (cavitation Index) 

H dimensionless film thickness, h/AR 

h film thickness, m 

L/D length to diameter ratio 

M bearing torque, N-m 

ifi x ,rtt 2 lineal mass flux, kg/m-s 

Pl power loss, N-m/s 

p fluid pressure, N/m? 

p a ambient pressure, N/m 2 

p c cavitation pressure, N/m 2 

p 0 dimensionless pressure, p(AR) 2 /(R 2 pw) 

R radius of shaft, m 

AR radial clearance, m 

t time, s 

U sum of the surface velocities In x-dlrectlon, m/s 

sum of the surface velocity vectors, m/s 
W load capacity, N 

w squeeze-velocity, m/s 

x coordinate along circumference, m 

Ax Incremental spacing along circumference, m 

y coordinate normal to x, z-plane, m 

z axial coordinate, m 
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az axial Incremental spacing, m 

0 liquid bulk moaulus, N/m 2 

T angular position of minimum film, rad 

e eccentricity ratio, e/AR 

ti dynamic viscosity, N-s/m 2 

p c fluid density within cavltated zone, kg/m 3 

<p angular coordinate relative to minimum film line, rad 

angular location of upstream cavitation boundary, rad 

<p e angular location of downstream cavitation boundary, rad 

e fractional film content In cavitation zone; density ratio (p/? c ) 

In full film zone 

ud orbital angular velocity of journal center about fixed point relative 

to housing center, rad/s 

u> s angular velocity of journal about Its own center, rad/s 

< 1 > attitude angle, tan (F^/Fr), rad 

t x shear stress for an oil element, N/m 2 

INTRODUCTION 

Machine elements In relative motion separated by a lubricating fluid can 
often be subjected to conditions which causes the fluid to cavitate. This 
cavitation can either be a result of (1) dissolved gas coming out of solution 
or (2) evaporation (flashing) of the fluid. Both types of cavitation are 
commonly observed In Journal bearings, squeeze-film dampers, connecting rod 
bearings, and pistons. It has a pronounced effect on their operation. The 
occurrence of cavitation In Journal bearings Is shown (1) to result In reduced 
power loss, friction coefficient, bearing torque, and load capacity. Dowson 
and Taylor (2), In an excellent review of cavitation, point out that 
cavitation need not have a deleterious effect on the load-carrying capacity of 
bearings. Horsnell's predictions (3) on load capacity for steady-state 
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conditions show that consideration of cavitation can lead to a threefold 
Increase over the predicted load If It Is not considered. 

In recent years, Imposing higher loads and speeds together with more 
complicated loading cycles resulted In bearings experiencing "cavitation 
erosion-damage" ( 4 ). This damage occurs under conditions of vapor 
cavitation. The vapor pressure for most lubricating oils Is very nearly 
zero. Very large Implosive forces are created and confined to a small area as 
the vapor bubble collapses. The net effect Is that metal Is hammered out of 
the bearing surfaces by a fatigue process, resulting In erosion damage to the 
bearing surfaces. This phenomenon arises under conditions of dynamic loading 
and Is frequently observed In main or crankshaft bearings In compression- 
ignition engines. 

Aside from the cavitation damage, there Is the damage that can occur due 
to ielf-exclted Instabilities that are encountered under dynamic loading. 

Those Instabilities can he manifested as a whirling or whipping motion (5) of 
the Journal center. Larg^ vibrational amplitudes can result In large forces 
being transmitted to the system. A designer can circumvent this problem, or 
at least the severity of It, by knowing the speed threshold at which it 
occurs. This Information Is often obtained by referring to "stability maps" 
(6, 7). Stability maps require the determination of hydrodynamic force terms 
coupled to the dynamical equations of motion. 

The calculation of the hydrodynamic force components Is dependent on the 
f^im model used, especially at high eccentricities. Many of the film models 
that are used to generate stability maps Involve theories that are overly 
restrictive (l.e., narrow bearing theory. Infinite bearing theory) and treat 
the cavitation In a very superficial way. Notable of these Is the w-fllm 
cavitation theory In which positive pressures extend through half the 
circumference of the bearing and the other half Is regarded as cavltated. 
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Specifying zero-pressure gradient and cavitation pressure at the boundary 
(l.e., Swlft-Stleber boundary conditions) represents the conditions for film 
rupture In a reasonable way but does not properly represent the conditions 
when the film Is reestablished. In actual practice, the Swlft-Stleber 
conditions are often compromised as a matter of convenience by neglecting the 
Implementation of the zero-pressure gradient condition. This Is much like the 
Gumbel conditions In which solutions wure obtained to the Reynolds equation 
and cavitation was determined by disallowing the existence of subambient 
pressures. However, dynamic loading causes changes In the local film 
thickness which leads to nonstationary cavitation - that Is, the cavitation 
boundary Is In motion. This motion Is manifested as growth anti collapse of 
the bubble as well as downstream transport from the minimum film position. 

The appropriate boundary condition may require a condition other than a 
zero-pressure gradient at the boundary. The studies of Olsson (8,9) suggest 
that the usual Swlft-Stleber film-rupture condition Is adequate for dynamic 
situations if the cavitation boundary moves at a speed that Is less than half 
the Journal surface speed. It seems highly unlikely that this condition would 
be met during the Initial stages of growth and the last stages of bubble 
collapse. 

Jakobsson-Floberg (X) and Olsson (8) formulated boundary conditions for a 
moving boundary that conserved mass within the cavltated region as well as at 
the boundary (commonly referred to as the JFO cavitation theory). It was 
assumed that liquid was convected through the cavltated region In the form of 
strlatlons extending to both surfaces In the film gap. There Is evidence that 
some liquid Is transported In the form of an adhered layer to the faster 
surface (10, H). Pan (II) has broadened the JFO theory to accommodate those 
situations In which the mass transport through the cavltated region Is not 
necessarily via liquid strlatlons. For moderately to heavily loaded bearings 
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(1.e. t the load carrying capacity » surface tension forces of the 
lubricant), the adhered layer can be neglected (lfi). The striated flow In the 
JFO theory Is necessarily a Couette flow because of a constant pressure 
assumption within the cavltatod region. Floberg's experimental findings (1) 
support the constant pressure assumption. However, Etslon and Ludwig (lj.) 
have measured pressure variations of the order of SO kPa inside a gaseous 
cavitation bubble. To resolve this apparent contradiction. It seems that one 
must make the distinction of whether or not the cavitation bubble Is gaseous 
or vaporous. Etslon (1J.) provides a plausible mechanism based on liberation 
and reabsorption rates between a gas and a Mquld. The same arguments for 
pressure variation within a vapor bubble would not apply (cf. H). In view of 
this discussion, the constant pressure assumption for vapor cavitation seems 
reasonable. Despite some of the controversy and certain lack of understanding, 
the JFO theory perhaps represents one of the best accounts of a dynamical 
theory to date for moderately to heavily loaded Journal bearings and/or 
dampers, 't Is an Improvement over the Swlft-Stleber conditions, even for 
steady-state solutions, because It provides for film reformation or flllback 
as well as film rupture. Both rupture and flllback require a knowledge of the 
pressure gradient and fractional film content at the Interface to determine 
Its location. Unlike a rupture boundary, the flllback boundary Is subjected 
to a pressure flow (l.e., nonzero pressure gradient). Furthermore, the 
fractional film content at the boundary Is determined by a residual fluid 
within the cavltated region that has been released at a rupture boundary 
(earlier In time) and governed by the fluid transport law. Which condition 
(flllback or rupture) prevails depends on the relative motion of the boundary 
with the motion of the convected fluid normal to the boundary. Typically, 
under steady-state cavltatlng conditions, the rupture boundary occurs along 
the upstream boundary and the flllback boundary occurs along the downstream 
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end. Under dynamic loading and nonstationary cavitation, rupture or flllback 
conditions may be required at either the upstream or the downstream boundary. 
To effect these conditions In a computational algorithm, the programming task 
Is exceedingly tedious, and this discourages their implementation In current 
practice. 

Elrod and Adams (14) Introduced a computational scheme that mimics the JFO 
theory. It avoids the complex programming required to trace the moving 
boundary between grid points and evaluating the required pressure derivative. 
Later, Elrod (IJ.) modified It and presented It In more detail. The algorithm 
Incorporates a switching function that "switches out" any pressure flow within 
the cavltated region. This automatically Introduces a cavity conforming to 
the requirements of mass continuity and the JFO theory. Elrod has compared 
results generated by the algorithm with experimental results of Lundholm ( Tj6) 
for steady-state operation of a circumferentially fed journal bearing. The 
results (l.e., cavitation extent, load capacity, and attitude-angle versus 
eccentricity) were quite good. To the author's knowledge, a limited few 
(12,18) have used the algorithm for steady- state applications. No one to date 
has explored the full Intended capabilities of the algorithm, that Is, a 
dynamically loaded application. 

It is the purpose of this pap&r to analyze a dynamically loaded Journal 
bearing through Its entire periodic motion. The prescribed motion of the 
Journal relative to the bearing will be undergoing circular whirl to conform 
with the experimental work of Jacobson and Hamrock (13,19). Comparison of 
predicted bubble life to that measured by Jacobson and Hamrcck (13) will be 
discussed. Furthermore, results obtained from using an often used 
nonconservative theory (l.e., pseudo-Gumbel BC) are compared with the Elrod 
algorithm. 
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BACKGROUND THEORY 


Elrod and Adams (14) Implemented a moving boundary scheme that avoids 
Interface complications. The following Is a review of some of the underlying 
factors that lead up to the algorithm which was later modified by Elrod (1_5) 
and used here. The conservation of mass can be written as 

♦ 3 • m « 0 (1) 


where m represents the lineal mass flux and Is given by 


m * p 


h? 

2 


~ 12 V vp 


( 2 ) 


Substituting Eqs. (2) Into (1) leads to the Reynold's lubrication equation. 

This equation has been made applicable to the cavitation region as well as the 
full film region by Incorporating a switch function to automatically satisfy 
the boundary conditions at a moving interface. Furthermore, the fractional 
film content e has been made the dependent variable. This required giving 
e a dual Interpretation. That Is, In the full film region, 6 represents 
the mass content of the film that exceeds the content that would exist if the 
pressure were at cavitation pressure p . In other words, 

e - p/p c (3) 

where p Is the density of the liquid at the pressure p . Furthermore, 
c c 

p and p are related through the equation for the liquid bulk modulus 
according to 


or 



p * P c + 0 In 8 


8 


(4) 


e > l .0 


(5) 



In the cavitation region (s < 1 . 0 ) © determines the mass content (;> oh) 
which can manifest Itself In the form of a smeared mass or striated flow 
extending to both surfaces In the film gap. The fluid transported through the 
cavitation region In the form t:' an adhered film can be neglected for heavily 
loaded conditions. 

A universal differential equation Is made possible by linking the 
solutions of the full film region with the solutions In the cavltated region 
via a single dependent variable ft. However, a cavitation Index, or switch 
function g, was Included so that the resulting POE would be consistent with 
the uniform pressure assumption within the cavltated region. The switch 
function Is defined from a knowledge of ©. Thus, 


g « 1 ©>1.0 


and g Is mad* a factor of the pressure gradient term In Eq. (2) so that the 
flow Is strictly Couette In the cavltated region. Expressing the lineal mass 


flux m In terms of © and g and substituting It Into Eq. (1) result In the 
universal differential equation obtained In Ref. 14; l.e.. 


♦ | • v(©h) - d . ^ Bg(e)ve (7) 

In the full film region the solutions of e together with Eq. (5) 
determine the pressures. In the cavitation region, g » 0, and Eq. (7) becomes 

| • $(©h) - 0 (8) 

which governs the transport of the fluid through the cavitation region. 
NUMERICAL ALGORITHM 

The numerical scheme Is based on the Elrod algorithm (15) In which mass Is 
conserved throughout the computational domain. This method combines a control 
volume (CV) approach to deriving the finite-difference equations with the use 
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of a switching function to automatically satisfy the boundary conditions at a 
moving Interface, a CV Is constructed about each nodal point In space 
(Fig. 1) to Implement the conservation property In the differencing scheme. A 
computational cell of this sort requires evaluations of variables and their 
derivatives to be carried out at the cell faces. Variables at the cell faces 
are determined to be the simple arithmetic mean between Its adjacent nodal 
values at a time t. The derivatives are evaluated by centered differences. 
With this preamble, the conservation law Is stated as follows: 


The total Increase of 
mass (prhe) In the CV 
over a time At 


The net lineal mass flux 
of (p c he) Into the CV by 
convection plus the net 
lineal mass flux of (p c he) 
Into the CV by diffusion 


This translates to 


*(/> c eh) 

(A X AZ) — g* - Am^AZ ♦ Am^Ax (9) 

which Is equivalent to Eq. (1). The treatment of the lineal mass flux Is the 
crux of the Elrod algorithm and the essence of It Is as follows. Considering 
for the moment one direction, namely x, then 


m x “ convection * ^"Vpressure 
Convective Contribution 

The switch function g Is judiciously applied to the convective term at 
the upstream cell face to give 

"x-Ax/2 * 2 [ 9 -l h -l (1 - 9 -l> * «-l h l * V-l (h o - h -l )/2 ] <’ 

The subscripts refer to the nodal position relative to the CV (see Fig. 1). 
Eq. (11) Is well defined In that It Is self-consistent with the CV approach 
and retains the physics as well. The lineal mass flux at the downstream face 
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Is obtained by suitably reassigning the subscripts of Eg. (11). Thus the 
Incremental change to the lineal mass flux due to convection Is 

< 4 Vco„v *[¥ h -l (1 - - V' - V 6 . * !s f sl <* - V 

0 A h_ g.h.9- I 

♦ -f - 2 <s_, -?♦«,)- “V^J < 12 > 

An attractive feature of this expression Is that In the full film region (all 
g > 1) the factor d(hd)/dx that appears In the convective term Is central 
differenced, retaining second-order accuracy. In the cavitation region, the 
algorithm accounts for the mass transport consistent with conservation 
properties by Including e In the derivative (l.e., d(he)/dx). The role 
of the switch function automatically effects an upwind differencing scheme for 
the evaluation of this term, thereby retaining properly posed conditions at 
the boundaries as time Is advanced. 

Pressure Gradient Contribution 

The Polseullle or pressure gradient contribution to the lineal mass flux 
at the upstream face Is 


LW f v 

<V*x/2>pres, “ ' ■<•>(#) <13) 

where the bar represents the average value at the Indicated cell race. The 
pressure Is determined from e by the approximation to Eq. (5), 

p a p c ♦ 0(6 - 1) (14) 

Central differencing the pressure at the upstream cell face and making use of 
the previous equation one obtains 


(m x-A/2Vess 


h> »« - 9 o< 9 o - 3 > - «.i < 9 -i - n 

12 p 0 AX 


(15) 
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The lineal mass flux due to diffusion that Is leaving the CV at the 


downstream end Is 


^ m x*Ax/2^press 


h 3 
h p 


12S 0 


g l (e l 


- 1) - 9, 


AX 


(8 o - 


1) 


( 16 ) 


In this analysis we do not need to represent h by an averaging method. 
We can evaluate It exactly. Thus, the net lineal flux Into the CV by 
diffusion Is 


(AfO 

x'press 


r 3 

pjJB hf 

12p l 


1/2 9 -l ( 8 -1 


~ ] > - (h ?i/2 * h ?i/2 )9 o (8 o ~ 1} 


Ax 


- ’’l 

AX J 


(17) 


If all three points are within the cavitation zone, then Am x - 0. This Is 
consistent with the zero pressure gradient assumption In that region. If all 
three points are In the full film zone, then Eq. (17) reduces to a central 
differencing scheme for the pressure gradient. There are numerous 
Interpolatlve combinations that occur at the boundary using Eqs. (12) and 
(17), all of which make up the boundary conditions automatically and 
consistent with mass conservation from grid point to grid point. 

Time March 

The time march related to Eq. (1) must now be considered. In this study, 
an alternating direction Implicit (AOI) scheme was applied. The Implicit 
Euler method was used to advance the time at each half time step since this Is 
known to be unconditionally stable (20) Insofar as the ordinary differential 
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equation Is stable. For the first half time step, a differencing along the 
circumference Is performed; that Is, 


(he)*’ 4 *' 2 - (he)* . fj- (- 

c 


• t+6t/2 
*■* 


•?*) 


or 


( 18 ) 


ihe - 


2p. AX/ 


1 \t*At/2 


he 


2, c az ) 


During the next half time step, the axial direction Is differenced giving 


/ ’ UAt/2 * t*At' 

(he) . - (he) • 2T \IT * ST , 


or 


( 19 ) 


he - 


Am \t»At 


2 P r AX/ 


(■• • ft *) 


: \t+At/2 


Solutions 

The previous equations lead to two systems f equations that can be 
represented by a periodic trldlagonal matrix B p In the circumferential 
sweep and a trldlagonal matrix B for the axial sweep. In other words, we 
can represent Eq. (18) as 


2VJ e*’ 4 *' 2 . c J> * 

“ P 1 


( 20 ) 




Because of the periodic or wrap-around boundary conditions, B p J takes the form 

I b, c, a, I 


.'3 


2 b 2 c 2 0 
a 3 b 3 c 3 


0 


s b k 


( 21 ) 
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The solutions of 6 along the circumference are found by performing a 

Gauss-Jordan reduction on B using a maximum-pivot strategy to reduce the 

P 

error. The solutions obtained by differencing along the axial direction are 
more readily obtained by using a nonpl voted Gaussian elimination procedure 
(trldlagonal solver). Since the journal and housing are aligned, the 
pressures must be symmetric about the axial center. Consequently, the 
calculations are only made over half the axial length of the housing. 

Switch Function 

The switch function was updated after each half-time step. Occasionally, 
during the collapse of the vapor bubble, liquid was piling up at the 
boundary. It appears that the switch function was not accommodating the 
movement of the boundary adequately for that time step (£ 1 ). Consequently, 
the switch function was Immediately updated and the calculation was reiterated. 
RESULTS AND DISCUSSION 

A numerical scheme based on the Elrod algorithm was devised to study vapor 
cavitation In a submerged Journal bearing for both steady and dynamic 
loading. The method conserves mass throughout the computational domain 
Including the region of cavitation. In this Investigation, the algorithm was 
checked against published experimental data. Furthermore, an analysis of a 
journal whirling In a circular path through one orbit was then performed. The 
consequences of applying a nonconservative theory (pseudo-Gumbel boundary 
conditions) In lieu of the conservative Elrod algorithm are determined. 
Experimental Comparison 

Steady state . - lakobsson and Floberg (1) have measured pressures within a 
submerged Journal bearing operating at steady-state conditions with 
cavitation. The operating conditions for the particular experiment chosen as 
a basis of comparison are given In Table I. These conditions were used as 
Input to the computer code. Figure 2 Illustrates the bearing configuration 
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along with the predicted steady-state pressure distribution that resulted from 
the author's numerical calculations using the Elrod algorithm. The Journal Is 
rotating In a counterclockwise direction. Thus, the entrainment velocity 
through the minimum film Is In a counterclockwise direction, and by convention 
considered positive. The line of centers forms a horizontal which extends 
through the position of minimum film thickness. The positive pressures seen 
In the three-dimensional plot are generated In the converging clearance region 
below this horizontal line. In the direction of rotation, the pressures 
become subambient In the diverging clearance region above the horizontal. The 
upstream and downstream meniscus of the cavitation boundary Is outlined 
there. The corresponding extent of cavitation Is shown by the outline of 
asterisks In the three-dimensional plot. 

Pressures were measured along the circumference of the bearing at two 
axial positions relative to the axial center In Ref. 1. Figure 3 compares the 
predicted pressures using the Elrod algorithm with these experimental data. 

The vertical axis Is a measure of the dimensionless pressure p Q and the 
horizontal axis represents the circumferential position (In rad) relative to 
the minimum film line. In general, the agreement appears quite good. The 
extent of cavitation determined experimentally can be seen to be greater than 
the numerical solution at both positions. Furthermore, there appears to be a 
slight discrepancy at the position of peak pressure. These discrepancies did 
not appear to adversely affect the agreement between the predicted and 
experimentally determined load capacity. The predicted load (2147 N) came to 
within 4.6 percent of that measured In Ref. 1 (2250 N). 

Dynamic loading . - For practical considerations, very little work has been 
done on submerged journal bearings In dynamically loaded conditions. This Is 
because submerged journal bearings with ambient supply pressure require the 
existence of subambient pressures to remain operative. Otherwise, the flow 
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balance Into and out of the bearing cannot be maintained, especially If a 
portion of the duty cycle Is steady state. The fluid mi'st be degassed so that 
the gaseous cavitation that might otherwise form at near ambient pressure does 
not preclude the formation of subambient pressure. 

Jacobson and Hamrock (13,19) have studied vapor cavitation In a submerged 
journal bearing. In their studies, a high-speed motion picture camera was 
used to capture the cavitation formation and collapse under dynamically loaded 
conditions. Figure 4 Illustrates the configuration used In the tests. The 
bearing (housing) vibrated In a circular motion keeping the axis of vibration 
parallel to the axis of rotation. The journal could be rotated about Its own 
fixed center In space. The dynamic conditions chosen as a basis of comparison 
are stated In Table I. The Important result of Ref. 13 Insofar as this study 
Is concerned Is the observance of the complete life cycle of vapor 
cavitation. High-speed photography revealed that, for the conditions listed 
In Table I, the vapor bubble remained visible for at least 23 ms. This 
represented approximately one-third of the full dynamic cycle of the bearing. 

The predicted life of the vapor bubble from the computer code was 32 ms. 

A part of this difference arises because this measure of bubble life 

Inherently contains a certain time duration for which the computer Indicates 

cavitation although In reality It would be Invisible to the eye. This Is 

probably not the major contributing factor to the discrepancy, however. The 

7 2 

pressures calculated for these conditions were of the order of 10 N/m . 

A plastic (polymethylmethacrylate, PMMA) was used as a housing material for 
viewing purposes. It Is thus most likely that deformation effects are 
responsible for the discrepancy not accounted for by the numerical method. 
Furthermore, the experiment allowed the eccentricity to reach a maximum value 
of 1.0, whereas a certain amount of numerical Instability was observed when 
the eccentricity exceeded 0.98. Consequently, the Input conditions for the 
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computations were not Identical to the Input conditions of the experiment. 
Future experiments are planned so that comparison between theory and 
experiment can be more conclusive. 

Comparison of Two Theoretical Models 

Motion . - The conservative cell method (Elrod algorithm) 1$ compared with 
a nonconservative scheme (psuedo-Giimbel BC) to analyze the effects of vapor 
cavitation In a dynamically loaded Journal bearing. The prescribed motion 
(Table I) of the Journal center was In a clockwise circular orbit (-92.7 rad/s) 
about a point fixed In space relative to the bearing center. The Journal spun 
(-39.5 rad/s) In the clockwise direction about Its own axis, which was 
considered parallel to the axis of the bearing housing. The net result for 
this motion Is a positive entrainment velocity (Fig. 4) which, by convention. 

Is In the counterclockwise direction. A complete orbit took a total of 
66.7 ms. The Instantaneous eccentricity < Increased from an Initial value 
of 0.1 to a maximum of 0.8 during the squeeze action of the Journal (first 
half cycle). Ourlng Journal separation (remaining half cycle), it returned to 
Its Initial value of 0.1 . 

Cavltlatlon and pressure distribution . - The results obtained by using the 
Elrod algorithm (Fig. 5) are compared with those obtained from the often used 
psuedo-Gumbel (see the appendix) boundary conditions (Fig. 6). Figures 5 and 
6 Illustrate the similarities and differences of the pressure distribution as 
they occur throughout the entire orbit. Figures 5(a) and 6(a) represent the 
position of the Journal within the housing and the associated pressure 
distribution at the Initial Instant In t A me. The pressure buildup due to the 
combined squeezing and sliding motion of the Journal Is shown for the first 
half cycle In Figs. 5(a) to (d) and 6(a) to (d). It should be noted when 
comparing parts (c) and (d) that the maximum peak pressure does not occur when 
the eccentricity Is a maximum as one might first expect. This can be 
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explained by referring to Fig. 7 and noting that the squeeze velocity 

approaches zero and the sliding velocity approaches a minimum at the end of 

the half cycle. This offsets the effectiveness of generating Increasing 

pressures due to a decreasing film. The onset of cavitation occurs between 

parts (c) and (d), 27 ms Into the cycle. The squeeze component Is known to 

suppress the conditions for cavitation (22). Note that In Figs. 5(d) and 

6(d), the pressure distribution and the extent and position of the vapor 

bubble are distinctly different. The conservation algorithm of Elrod predicts 

a sightly higher pressure peak and a smaller vapor bubble, which begins about 

15° further downstream from the minimum film line. Proceeding from (d) to 

(e) shows the Journal /near the minimum film line) separating from the 

bearing. The Initial stages of separation creates a suction effect, causing 

the pressure hump to dissipate and the vapor bubble to expand. Note that 

parts (b) and (e) were both generated at the same value of eccentricity. The 

essential difference Is that the squeeze velocity has a different sign; that 

Is, when the squeeze velocity Is negative (the clearance Is decreasing), the 

motion of the Journal Is producing a squeeze effect (part (b)). In part (e), 

the squeeze velocity Is positive (l.e., the clearance Is Increasing) and 

produces a suction. In comparing Figs. 5(e) with 6(e) the difference In 

bubble size has become even more noticeable. Here we see In both cases that 

the vapor bubble has actually crossed the minimum film line and been drawn 

Into the converging clearance space of the bearing. This effect Is greater 

for the nonconservative theory. As the Journal continues to pull away from 

the bearing, the Increased clearance means that a much greater Polseullle side 

« 

flow Is present to cause the collapse of the bubble. This Is because the 

3 

Polseullle side flow Is proportional to h . In both theories, the bubble 
drifts downstream from the minimum film line and collapses. The Elrod 
algorithm predicts a much longer bubble life (8.4 ms) with an oblong shaped 
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bubble that actually crosses the maximum film line before collapsing. This 
collapse occurs about 7 ms Into the succeeding orbit. 

Bearing film calculations . - The effects of boundary conditions on bearing 
film calculations are shown In Fig. 8. The calculations derived from the 
nonconservative model (dotted lines) are superimposed on the curve derived 
from the conservative model (solid line). The load capacity and load 
components are shown In Figs. 8(a) to (f) as a function of the eccentricity 
during a full period of shaft whirl. The letter callouts on the plots can be 
used as chronological Indicators since they refer to the parts In Figs. 5 and 
6 chronologically. The load-eccentricity plots might well be described as 
"Journal bearing hysteresis loops." For Increasing c the W-« curve Is 
entirely different from that for decreasing «. That Is, proceeding from 
part (a) to part (d), c Is Increasing and the squeeze velocity Is always 
negative. Ourlng this Initial half-cycle, cavitation Is for the most part 
absent. For decreasing «, the squeeze velocity Is always positive and 
cavitation Is almost always present. The drastic loss In load capacity In 
going from part (c) to (d) Is a result of the Journal motion (Fig. 7) and the 
development of cavitation close to the minimum film line. The negative 
pressures In the diverging clearance tend to cancel the positive pressures In 
the converging clearance. However, the motion Is perhaps the major factor 
Influencing the drastic load loss. 

As should be expected, differences from the boundary conditions appear 
only when cavitation Is present. At the higher eccentricities 
(0.5 < c < 0.75), the largest difference In load capacity Is about 
10 percent. The radial load component accounts for the largest part of this 
deviation. Again this Is because the vapor bubble Is located close to the 
minimum film. In fact. In the discussion of Figs. 5 and 6 It was pointed out 
that the bubble actually penetrated the region of converging film. Note that 
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a larger bubble leads to a reduction In the load components. The reverse 
effect Is noted for the total load capacity. This Is because the magnitudes 
of the load components are the Important factors and not the relative values 
when determining total load (see the appendix). As the eccentricity continues 
to decrease, the dotted and solid lines tend to converge, except In the case 
of the tangential load. In this plot, these two lines continue diverging 
until the journal begins Its next orbit. In the nonconservative theory the 
bubble has collapsed just 1 ms before the journal returns to Its Initial 
position. By contrast, the vapor bubble plays a very active role In the load 
characteristics well Into the next orbit, according to the conservative 
theory. As the bubble penetrates the maximum film, the radial load once again 
becomes the dominant factor In the deviations due to boundary conditions. 

These deviations represent as much as a 20-percent difference In the total 
load capacity. The dramatic dip In load during the Initial stage of the 
second orbit and each succeeding orbit Is due to the sudden collapse of the 
vapor bubble. The hysteresis loops are repeatable Insofar as the circular 
orbit of the journal Is maintained. 

The attitude angle Is an angular measure between the load line and the 
line of centers. It Is defined In the appendix along with the film force 
calculations. The effect of the boundary conditions on the attitude angle Is 
seen only when the eccentricity Is decreasing. This effect results from the 
differences In extent and location of the bubble. 

The bearing torque and power loss In the bearing show no boundary 
condition effect except perhaps on a fine scale. The only differences In 
making these calculations are In the handling of the shear stress In the 
cavltated region. Since pressure 'gradients are assumed not to exist there, 
the shear stress arising from the Couette flow must reflect any possible 
differences. This calculation Is directly proportional to the fractional film 
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content In the vapor region and the extent of the vapor region. In the case 
of the nonconservative theory, the fractional film content 1$ not calculated. 
Therefore, It must either be assumed to be entirely vapor or filled with oil. 
These calculations were made assuming the zero pressure region, which has a 
linear velocity distribution. Is filled with oil. According to Figs. 8(e) and 
(f), this contribution was negligible for the scale shown. 

CONCLUDING REMARKS 

A theoretical Investigation Is made of the evolution of a vapor bubble for 
a submerged journal bearing under dynamically loaded conditions. The solution 
to the Reynolds equation Is determined numerically using a control volume 
method (Elrrd algorithm). This method conserves mass throughout the 
computatloral domain Including the liquid-vapor Interface which may or may not 
be In motion relative to the minimum film line. An AOI method Is used to 
effect the time march. 

Excellent agreement was found with the experimental work of Jacobson and 
Floberg (1) for stationary cavitation. A prediction of 32-ms bubble life for 
nonstationary cavitation was obtained from the dynamic theory for conditions 
(13) listed In Table I. Using high-speed photography for the same set of 
conditions, Jacobson and Hamrock observed a bubble life of 23 ms. This 
discrepancy was attributed to deformation effects as well as some uncertainty 
In correlating actual life with duration of visibility. Furthermore, the 
maximum eccentricity In the experiments was 1.0. Numerclal Instabilities 
prohibited eccentricities larger than 0.98. Future experiments are planned to 
preclude some of these problems so that comparison between theory and 
experiments can be more conclusive. 

A comparison study was performed to determine some of the consequences of 
applying a nonconservative theory to a dynamic problem. A complete dynamic 
cycle of a journal whirling In a circular path was chosen for the basis of 
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comparison. Significant differences were observed In the load components near 
the end of the cycle. As much as a 20-percent difference In load capacity was 
observed. Differences In attitude angle were brought about by differences In 
extent and location of the bubble. Bearing torque and power loss were not 
affected by applying the pseudo-GUmbel boundary conditions In lieu of the 
dynamic theory. 

In conclusion, the consequences of applying a nonconservative theory are 
noticeable only In the load components at the tall end and beginning of the 
Journal orbit. The ease with which the Elrod algorithm can be Implemented Is 
an attractive feature. However, the present code consumes roughly 2 to 3 
times the computational time of the GUmbel solutions when using a Gauss-Seldel 
Iterative SOR scheme. Faster codes are available (l.e., fast Fourier 
transforms and multigrid codes) that could conceivably reduce the 
computational time by an order of magnitude. This certainly would make the 
algorithm a viable film model for rotor dynamic and squeeze-film damper 
applications. 



APPENDIX 


Film Force Calculations 

The film forces are all determined here In much the same way as they were 
derived In Ref. 1. The tangential load component Is a force acting normal to 
the line of centers at the minimum film position: 

F ■ - J/ S p sin (p R dtpdz (Al) 

The radial load acts along the line of centers: 

Fr - JJ S P cos <9 R d<pdz (A2) 

The total load and attitude angle can then be determined: 


w .y 


f r f F? 
R <p 


♦ - tan-V^) 

The friction force is determined from the shear stress t x : 


% ■ [L ’xl R d «' d2 

JJ* v»h 


y-h 

For a film thickness h the shear stress Is 


h io ye 
x|y*h " 2 ax * u *- 

where . . 

h a AR (1 - e COS <p) 

H * 1 - e COS <p 
If one Integrates the pressure gradient term by parts, then 

\ ■ f/Hl; di *Jj *§*■ p s,n * d,pdz * “m fj * d * dz 

dz - w\'ff 


AR 

2 


AR H R d<pdZ 


(A3) 

(A4) 


(A5) 


(A6) 


(A7) 
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All Integrations ware performed using a Simpson's rule for single as well as 
double Integrations. 

Pow er loss 

The power loss can b e obtained from the shear force, 1. e. f 

P, • F U (A8) 

. B ear in g Tprqg e 

The bearing torque Is given as 

H.FR (A9) 

v 

Pseudo-Gumbel Boundary Conditions 

The Gumbel boundary conditions are often used as a substitute for the more 
rigorous Swift-Stleber rupture boundary conditions for reasons of 
computational economy. Normally the Gumbel conditions are Implemented by 
first solving the Reynolds equation without any regard to cavitation. All 
subambient pressures are disallowed by setting them equal to the ambient 
pressure. This serves as the condition for cavitation. Note that the 
specification of the pressure gradient at the boundary Is neglected. In this 
study, the nonconservative theory was not strictly In accord with Gumbel; that 
Is, fubamblent pressures were allowed, and the vapor bubble was determined by 
disallowing pressure* less than the cavitation pressure. 

The numerical procedure used to Implement the Gumbel conditions was the 
commonly used Gauss-Seldel Iterative technique. This was done to provide a 
point of reference In assessing the penalty (If any) of Invoking a more 
accurate theory. 
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TABLE I. - OPERATING CONDITIONS 


Steady-state 

, Dynamic 

Reference 1 

Reference 13 

Present study 

1.455xlO' 4 m 

5.0xl0 -4 m 

5.0xl0 -4 m 

0.050 m 

0.0425 m 

0.0425 m 

4/3 

1/4 

1 .0 

0.61 

0.32-1 .0 

0.1-0. 8 

-48.1 rad/s 

-19.5 rad/s 

-19.5 rad/s 

0.0 rad/s 

-92.7 rad/s 

-92.7 rad/s 

1 .72x1 0 9 N/m 2 

1.72xl0 9 N/m 2 

1.72X10 9 N/m 2 

0.0127 N-s/m 2 

0.066 N-s/m 

0.066 N-s/m 

0.0 N/m 2 

1.0133x10 N/m 

1.0133x10 N/m 

-72139.79 N/m 2 

0.0 N/m 2 

0.0 N/m 2 
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Figure 1. - Control volume CV at point x. 



Figure 2. - Pressure distribution and bearing configuration using experi 
mental conditions of Floberg (ref. 1). 



(a) Axial position from centerline (L/5D). (b) Axial position from centerline (3L/5D). 

Figure 3. - Comparison of experimental to predicted pressures (dimensionless) along circumference. 
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(f) Bubble collapse (t * 65.0 ms ; e * 0.141. (el Bubble expands (t ■ 53.4 ms ; e * 0.501. <d) Begin separation (t * 33.6 ms ; £ * 0.80. 

Figure 6. - Pressure distribution and bearing configuration for a full period of shaft whirl using pseudo-Gumbel boundary conditions. Figures (al to If) viewed in clockwise order are 
consecutive in time. 
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